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Abstract 

The magnetic properties of a system of coexisting localized spins and conduction electrons are 
investigated within an extended version of the one dimensional Kondo lattice model in which 
effects stemming from the electron-lattice and on-site Coulomb interactions are explicitly included. 
After bosonizing the conduction electrons, is it observed that intrinsic inhomogeneities with the 
statistical scaling properties of a Griffiths phase appear, and determine the spin structure of the 
localized impurities. The appearance of the inhomogeneities is enhanced by appropriate phonons 
and acts destructively on the spin ordering. The inhomogeneities appear on well defined length 
scales, can be compared to the formation of intrinsic mesoscopic metastable patterns which are 
found in two- fluid systems. 

PACS numbers: PACS No. 71.27.+a, 71.28.+d, 75.20.Hr 
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The interplay of spin, charge and lattice degrees of freedom has been investigated inten- 
sively in many transition metal oxides and especially in perovskite manganites, which have 
recently attracted new interest due to the discovery of colossal magnetoresistance (CMR). 
The initial understanding of the properties of manganites was based on the double-exchange 
(DE) mechanism within the Kondo lattice.— However, the new experimental findings^ have 
revealed that this approach is incomplete and has to be extended to account for effects 
stemming from the lattice in order to understand the doping dependent phase diagram and 
the richness of phases that appear. In the following we model these complex systems within 
the Kondo lattice model (KLM), admitting for ferro- and antiferromagnetic couplings, and 
including explicitly the interaction with the lattice degrees of freedom. 

The KLM considers the coupling between half-filled narrow band (localized d or /) 
and conduction electrons. Even though studied intensively for the last two decades, the 
understanding of the KLM remains incomplete. Only in one dimension have numerical 
simulations 3 and bosonization techniques^ have been carried through which admit predic- 
tions about the phase diagram of the KLM. No investigations exist for the case where the 
KLM is extended to account for contributions stemming from the phonons, which is of spe- 
cial relevance to CMR materials. In particular, the small doping regime of these systems, 
which are ferromagnetic at low temperatures, seems to be appropriate to be modeled within 
the KLM extended by interactions with the lattice. In the following we present bosonized 
solutions of the KLM where on-site Coulomb and phonon contributions are explicitly in- 
cluded. This "extended" KLM model allows spin-and magnetoelastic-polaron formation, 
which we believe are of major importance in understanding these complex materials. 

The Hamiltonian of the KLM in the presence of on-site Coulomb interaction reads: 

Hklm = - i E(i c j+V + hx -) 

+ J Y, s d,i-s c j + u n iA n i,i > (!) 

3 j 

where t > is the conduction electron hopping integral, Sdj = § J2a,a' ^dj^^cryc^j ^,, 
S c ,j = \ J2a,a' Cj^v'vyCj^' an d cr are the Pauli spin matrices. Fermi operators c^, 
with subscript d refer to localized ci-spins, while those not indexed refer to the conduction 
electrons. The on-site Coulomb repulsion is given by the Hubbard term proportional to 
U. In the CMR materials the localized states are represented by the threefold degenerate 
Mn t2 9 (i-electrons with total spin 3/2. However, for reasons of transparency, the localized 
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spin is approximated here by spin 1/2, since the properties of the model are qualitatively 
independent of the magnitude of the localized spins. In the following the the Kondo coupling 
J is measured in units of the hopping t and both cases, antiferromagnetic (J > 0) and 
ferromagnetic ( J < 0) couplings, will be considered. The conduction band filling is given by 
n = N c /N < 1, where N is the number of lattice sites and N c is the number of conduction 
electrons. To be able to understand the wide range of properties of the CMR materials, we 
also allow for the number of impurity spins, Nd, to vary, in such a way that Nd/N < 1. 

In principle, the electron-phonon coupling could be of either inter-site (Su-Schrieffer- 
Heeger (SSH)S) or on-site (Holstein 7 ) character. Since we found the SSH-coupling to be 
irrelevant to forward scattering processes, its influence will not be discussed in the following; 
only terms arising from the on-site couplings,— i.e., J2j a( lj n ji with coupling constant a and 
displacement qj will be included. The bare lattice Hamiltonian is: H\ att . = ^2j\p 2 J2M + 
Kq 2 /2, where pj is g^'s conjugate momenta, K the harmonic coupling and M the ionic 
mass. 

The underlying bosonization scheme follows standard procedures^ by first decomposing 
the on-site operators into Dirac fields, c x a ~ J2 T e%kFX ^*T, (?{%)■, where Izf = nn/2, with spinor 
components r = ± (+/- being the right/left movers) and kp = nn/2. Next we bosonize the 
Dirac fields with \I/ ri(T = exp(i$ TiCr ) / V27rA, where 1/A is the ultraviolet cutoff. For the scalar 
Bose fields, $ Ti(T (x), and their conjugate momenta, il r o -(x), $ T)(T (a;) = Jf^ da;TI T)0 .(V), are 
used in standard Mandelstam representation by means of which a momentum cutoff via the 
Fourier transform is introduced A(/c) = exp(— \\k\/2). If the distance between the impurity 
spins is larger than A, the electrons will behave as collective density fluctuations. 8 Thus, the 
Fermi fields can be represented in terms of density operators which satisfy Bose commutation 
relations: c TXa ~ exp {it kpx) expi{9 p (x) + rcf) p (x) + a[9 a (x) + T(p a (x)]}/2, where the Bose 
fields for v = p,a are defined by ip u (x) = i(tt/N) J2k^o e lfca; [^ + (A;) ± u_(k)]A( y k)/k, with + 
corresponding to the number fields ip v = <p v and — to the current fields ip v = Q v . The 
charge (holon) and spin (spinon) number fluctuations are defined as p T {k) = J2a Pr,a(k), and 
cr T (k) = Y.a a PT,a{k). All rapidly oscillating terms originating from e.g. backscattering and 
umklapp processes are neglected, since they contribute only at exactly half filling. 

The localized d electrons can neither be bosonized nor Jordan- Wigner transformed since 
no direct interaction exists between them. Using the continuum approximation for the 
phonon contribution, two components are relevant: a small momentum part $oG°) and a 
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rapidly oscillating term at 2k F , $ n (x), resulting from the splitting of the conduction band 
electrons into right and left movers. While the former contribution causes forward scattering 
and is best represented in momentum space, the latter one gives rise to backscattering and 
requires representation in real space. The transformed Hamiltonian thus becomes: 



If holes are present in the array of d-spins, all terms proportional to S are zero. The 
notations used in Eq. (J2J) are: the forward scattering Holstein electron-phonon coupling 
term H cl ~ ph = (a/VM)(V2/N)J2 P [p+(-p) + p_(-p)]$o(p); the bare lattice contribution 

m h = (i/2A0E P [ng(p) + ««(/>)] + hfdxK(*) + <4&Ax)] with ^ = ^ = Jk/m; 

and the standard spinon-holon term H eX = (vp/Air) Y,j v \J^t{j) + [dx(j>u(j)] 2 } with velocities 
Vp/a = v F [l ± U/ttv f =f a 2 /irKv F ] 1 / 2 . 

It is important to note that a renormalization of the spinon-holon velocities appears 
here due to the Hubbard and phonon terms which act oppositely on the corresponding 
velocities. While the Hubbard term leads to a localization of the spinons and an increased 
hopping of the holons, thus supporting a magnetic ground state, the phonons delocalize the 
spins, but localize the charges and act destructively on the magnetic properties. It is worth 
mentioning that the Hubbard term alone already suffices to establish two time scales for the 
holon-spinon dynamics, but an important renormalization of the critical properties of the 
system is achieved through the variable phonon coupling, which as will be shown below 
establishes the existence of a Griffiths phase. The competition between the Hubbard and 
the phonon term obviously vanishes for U = a 2 /K. 

In the following effects arising from the localized spin d impurities, double exchange (DE), 
the phonons and Hubbard interactions will be discussed in more detail. The localized spin 
d impurities act via double exchange (DE) on the hopping electrons so as to preserve their 
spin when moving through the lattice in order to screen the localized spins which are in 
excess of the conduction electrons, i.e. N > N c . This, in turn, leads to a tendency to align 
the localized spins and results in an additional screening energy for the conduction electrons. 

In order to gain a more transparent understanding of this complicated interplay, the 
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j^r £ {C«# ff 0')] + COS [ 2 ^' + Mm (<T*'U)S+ j + h.C.) 




(2) 
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model is investigated first for the case of two sites and one conduction electron,— next in a 
simple continuum approximation, and finally the full bosonized solution will be presented. 

In the case of ferromagnetic coupling (J < 0) the ground state energy is E 0i j <0 = 
-|J|/4 - t with wave function |^o)j<o = |^de)j<o,* = I ft»T*» fz °) + I ft* 0, ft*T*)> 
where ff z and | z refers to the z component of the impurity and conduction electron spins, 
respectively. Ferromagnetism arises here via an Ising type coupling, which allows for de- 
scription of the ground state within a simple semiclassical approximation.— For J > the 
situation is completely changed due to the singlet formation of local and conduction elec- 
tron spins. This causes a mixing of the total spin and an enhancement of the Hilbert 
space, where now 16 elements have to be considered. The ground state energy is given 
by £0, j>o = - J/4 - V J 2 + 2Jt + 4t 2 /2 with wave functions \ip )j>o oc \ipKs)z + [1/(^/4 - 
£o,j>o)]{| tzlzA* 0} + | t* 0, 1U*)-| tizUAz 0}-| 4*0, tstz}}, where the Kondo singlet 
\ipKs)z states are | f z [ z , f z , 0) - | f z 0) + 1 f z 0, f z i z ) - | f z 0, Jj z T z ). |^o)j>o involves 

six basis elements (the degeneracy is partially lifted by conduction electron hopping) and 
hence falls outside the four dimensional space needed to establish DE for J < 0. In order to 
invoke DE as well, all three spin directions, x, y, and z, have to be considered: \ipo)j>o oc 
[1 - 1/(J/4-£ ,j>o)] |Vks)*+[1/(J/4-£ 0i j>o)] {|V j de)j>o, :c + |V'de}j>o, 2/ + |V'de}j>o, z }, 

Where \4>De) J>0, a=x ov y = {I taia, fa 0) + | fta 0, ft a l a ) + | J| Q ta, 0) + | ^ 0, Jjata)}/^ 

and IV'de) j>o,a=z = I fzlz, fz 0} + I fTz 0, ft~«J, z ). In spite of this extra complication, it is 
apparent from the above that in both cases, J < and J > 0, spin polarons are formed. 

Going beyond the two site approximation, spin polaron formation can be derived di- 
rectly from the KLM Hamiltonian, which can be written as: H « Ei(cL c i+i a + h.c.) + 
J/^J2i( n i,] ~ n i,i)Sd,j- Here spin-flip interactions are neglected, since these require a much 
higher energy and are consequently unlikely to be of importance to our results. This sim- 
plified Hamiltonian, as compared to Hamiltonian Eq. can be solved when the elec- 
tronic wave functions are treated within the continuum approximation and in the limit 
N c /N « 1, a case which is relevant to small doping concentrations in CMR materials. 
The electronic wave functions, %p a {x) ) satisfy a standard nonlinear Schrodinger equation: 
d 2 ip a {x) + (Jm e ij2)\%l) a (x)\ 2 il) a (x) = 2m e \Eip a {x) {m c \ being the bare electron mass) with 



These soliton solutions correspond to spin domain walls of finite size (kink-antikink pairs) 
and lead to a gain in electronic energy of — o for antiferromagnetic coupling, and of +a for 



soliton solutions ip a {x) oc e iX sech(xw Jm el /4). 
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the ferromagnetic case. Physically the solutions resemble the dressing of the electron by a 
finite range of parallel (antiparallel) local spins and consequently represent polaronic type 
objects. From the previous considerations it can also be concluded that, when including the 
interactions with the phonons, the tendency towards charge localization is enhanced and 
increases this polaronic effect. Since the lattice also experiences a renormalization due to 
the coupling to the electronic degrees of freedom, substantial ionic displacement patterns 
will develop and the formation of magnetoelastic polarons takes place. Similar results are 
obtained by decoupling electronic and phononic degrees of freedom through a homogeneous 
Lang-Firsov transformation, where the localization stems from band narrowing. The lo- 
calization width (polaron radius) is characterized by a length scale proportional to 1/ y/j . 
This new length scale differs from the free conduction electrons mean free path and gives 
rise to competing time scales: slow motion of the polaronic carriers and fast motion of the 
free electrons thus inferring dynamics of two types of particles and a close analogy to a two 
fluid scenario. Since the polarons are in general randomly distributed within the local spin 
array, these states can be viewed as intrinsic inhomogeneities involving spin fluctuations 
and short-range spin correlations. In addition these new slow dynamics will exhibit a peak 
in the spin structure factor at 2kp — tt instead of the simple 2kp RKKY signal. A similar 
observation has also been made^ using numerical approaches. 

In order to investigate rigorously the ordering of the local spins due to the formation 
of polarons, we first apply, an infinite (to avoid truncation errors) unitary transformation, 
S, to the bosonized Hamiltonian, Eq. (0). The most effective form of S is given by: S = 
i(J/27t)^vf/v^. J2j ^o-(j) S&j, which couples the conduction electron spins directly to the 
localized spins. Secondly, we explicitly take into account the Luttinger liquid character 
of the Bose fields, i.e., use their non-interacting expectation values such that the effective 
Hamiltonian for the local spins is derived as: 

#eff = -tV^£ j dk C os[k{ 3 - 3 ')]K\k)Sl 3 Sl r 
+ ^ Y.{™AKm + cos[2A*j]}S^ 

" ^£sin[^(i)W2M^ ■ (3) 

3 

Here K(j) stems from the unitary transformation and counts all the S^/s to the right of 
the site j and subtracts all those to the left of j: K(j) = (J/2vp) X^i(#f — S , |,_;)- This 
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term gives the crucial difference between the Kondo lattice and dilute Kondo lattice, as will 
be explained in the following. The most important term in Eq. (j3J) is the first one, which 
shows that a ferromagnetic coupling emerges even in the dilute Kondo lattice model. This 
coupling is non-negligible for > N c and j — j' < A and its strength will decrease with 
the distance between impurity spins. Thus, A represents the effective derealization length 
related to the spatial extent of the polarons, i.e., the effective range of DE. Thus DE will 
vanish if the distance between the impurity spins is larger than A. In general A will depend 
on J, N c or even N d , but we will use its low density value: A ~ \J2/ J. Consequently, we 
approximate it by its nearest neighbor form: X = ( J 2 v 2 a /2tt 2 vf) Jo°° dk cos kA 2 (k). 

For the Kondo lattice model, K(j) is vanishingly small as the number of <i-spins to the left 
and the right of a given site j is the same. The effective Hamiltonian can thus be replaced 
by a random transverse field Ising model:— H c s = S^jS^j +1 — J2j hjS^j, where the 

ferromagnetic coupling strictly vanishes if i > A. The random fields, hj, are generated by 
(1 + cos[2/cpj]) at large distances, where cos[2k F j] oscillates unsystematically with respect 
to the lattice. The large values cos[2kpj] « 1 which are responsible for spin flips, are then 
widely separated and are driven by a cosine distribution similarly to spin-glasses.— If we 
have a small concentration of holes in the array of localized spins, then - opposite to the 
previous case - K(j) is non-vanishing since the hole spins are no longer equally distributed 
to the left and the right of a given site. This yields K(j) ~ (— l)i (J/2vf), which gives rise 
to a staggered field and antiferromagnetic ordering. 

Since our main interest here is to explore the occurrence of ferromagnetism in the pres- 
ence of the Hubbard and phonon terms, we focus on the transition between the param- 
agnetic and the ferromagnetic phase. This is controlled by a critical coupling J crit = 
(7r/4) sin(7rn/2) {1 - U/[2-Ksin(7rn/2)} + a 2 /[2nK sm(irn/2)}} 1 / 2 . For values J < J crit . a 
paramagnetic state exists which is dominated by polaronic fluctuations. For J > J cr i t . fer- 
romagnetism appears. The transition between these phases is of order- disorder type with 
variable critical exponent 5 = J C rit./J- It can be seen that, in accordance with our previ- 
ous observation, the Hubbard term stabilizes the ferromagnetic phase, while the phonons 
counter this and tend to increase the polaronic regime. This paramagnetic polaron state 
can be viewed as a Griffiths phase, since the critical exponent is variable and the spin-spin 
correlation function is given by: (£/x) 5 / 6 e~( 3 / 2 )( 27r2:r /£) 1/3 e~ x ^, where £ ~ l/<5 2 is the cor- 
relation length. At finite temperature the susceptibility in this phase is proportional to 
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T 2 <5 -1 (lnT) 2 , whiig specific heat follows a T 2 ^ dependence. This regime can be viewed 
as a paramagnet with locally ordered ferromagnetic regions, again manifesting the analogy 
to a two fluid picture. 

In summary, we have derived an effective Hamiltonian from a one-dimensional Kondo 
lattice model extended to include effects stemming fom the lattice and in the presence of an 
onsite Hubbard term, which accounts for the conduction electron Coulomb repulsion. The 
results are: i) A ferromagnetic phase appears at intermediate | J\ due to forward scattering 
by delocalized conduction electrons, ii) Ferromagnetism is favoured by the Hubbard term, 
while it is suppressed by the electron-phonon coupling, iii) The paramagnetic phase is 
characterized by the coexistence of polaronic regimes with intrinsic ferromagnetic order and 
ordinary conduction electrons, iv) In the paramagnetic phase, two time scales compete with 
each other - reminiscent of a two-fluid model - and the variability of the critical exponents 
suggests the existence of a Griffiths phase. The results are related to the small-doping regime 
of CMR materials which are ferromagnets at low temperatures, since here the coupling to 
the phonons has been shown to dominate the paramagnetic-ferromagnetic phase transition. 

It is interesting to note the discrepancy between infinite dimensional calculations and 
the present one dimensional result. Many calculations to model CMR— have been made 
in dynamical mean-field theory, which is an infinite dimensional approximation and there- 
fore incapable of capturing spatial inhomogeneities. In the present work we approach the 
CMR materials via a one dimensional approximation, but with techniques able to describe 
fluctuations of short-range order. Our results show that strong intrinsic spatial inhomo- 
geneities of Griffiths type dominate the behaviour of the Kondo lattice. Consequently the 
inhomogeneities exhibit clear statistical scaling properties as a function of the proximity to 
a quantum (order-disorder) critical point. The phonons enhance the inhomogeneities, which 
in a good approximation behave as a supercritical (metastable) phase of a two fluid model. 

Even though various bosonization schemes have been used for the one-dimensional 
KLM,^ non of the previous approaches took into account phonons and the possibility of 
diluting the array of impurity spins. The inclusion of phonon degrees of freedom has been 
shown to be relevant in creating local magnetic inhomogeneities. It is important to mention 
that the properties of the system are driven by intrinsic inhomogeneities. This means that, 
in a renormalization group approach, the dimensionality should not matter-2 Thus, simi- 
lar behaviour is expected in realistic two- and three-dimensions, which clerly merit further 
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detailed study. 
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